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ABSTRACT 

We examine the colour-magnitude relation of ~ 5000 Brightest Cluster Galaxies 
(BCGs) in the SDSS. The colour-magnitude and colour-velocity dispersion relations 
of the BCGs are flatter in slope, by factors of two or more, than those of non-BCG 
early-type galaxies of similar luminosity (M r < —22.5), while their g — r colours at 
the half-light radius are ~ 0.01 magnitude redder. 

We investigate and compare radial colour gradients (which are usually negative) 
in these BCGs and a sample of 37000 early-type galaxies, using a gradient estimator 
based on the ratio of de Vaucouleurs effective radii in g and r passbands, ^777^ — 1. 



The mean colour gradient of BCGs is flatter (by 23%) than for other E/S0 galaxies of 
similar luminosity. 

The colour gradients in early-type galaxies are stronger at intermediate luminosity 
(M r ~ —22) than at low or the highest luminosities, and tend to decrease with in- 
creasing velocity dispersion (a). In non-BCG E/SOs, colour gradients increase for larger 
effective radii (up to 10-12 kpc), and are negatively correlated with 10 log a + M r , the 
mass density, and stellar age. However, gradients can be reduced or inverted (positive) 
for post-starburst galaxies at the youngest ages. In BCGs, these trends are absent and 
the mean colour gradient remains at a relatively low level (~ 0.08 by our measure) 
whatever the other properties of the galaxies. The redder half-light radius colours of 
the BCGs can be explained by their slightly greater ages combined with flatter radial 
colour gradients. 

We discuss possible explanations in terms of spheroidal galaxies initially having 
a colour gradient positively correlated with luminosity and positively correlated with 
large radius and/or low density. Subsequently, elliptical-elliptical 'dry' mergers pro- 
gressively reduce colour gradients, towards a low but non-zero value. This has occurred 
a greater number of times during the formation histories of the most massive E/SOs, 
and to by far the greatest degree in the BCGs. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters 
- galaxies: evolution 



1 INTRODUCTION 

In this paper we investigate colour-magnitude relations, and, 
to greater depth, the properties of radial colour gradients, 
for early-type galaxies and in particular Brightest Cluster 
Galaxies (BCGs). Both can provide information on these 
objects' formation. 

Early-type galaxies (elliptical and SO) with higher lumi- 
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nosities are redder in colour than those of lower luminosity, 
following a relatively tight and linear colour-magnitude re- 
lation (CMR), or 'red sequence'. Gallazzi et al. (2005) per- 
formed a detailed analysis of 2 x 10 5 galaxy spectra from the 
Sloan Digital Sky Survey (SDSS), using a set of 5 spectral 
indices (D4000, H/3, R SA + H 7 a, [Mg 2 Fe] and [MgFe]') to 
estimate mean stellar ages and metallicities. Gallazzi et al. 
(2006) and Jimenez et al. (2007), from spectral analysis of 
the early-type galaxies, conclude that the slope of the red- 
sequence CMR is primarily caused by an increase in mean 
stellar metallicity with galaxy mass or velocity dispersion. 
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Mergers may then affect the CMR by building more lumi- 
nous galaxies with the same composition as their less mas- 
sive progenitors, which would tend to flatten the CMR slope 
at L > L* luminosities (eg. Skelton, Bell and Somerville 
2009). 

Age also has a (secondary) role in explaining the colour 
trend - the most massive E/SOs formed early and rapidly (in 
< f Gyr) , while those of lower luminosity formed stars over 
longer timescales, giving slightly younger flux-weighted ages. 
There is also evidence that E/SO galaxy formation occurred 
slightly earlier in dense environments (eg. Sheth et al. 2006, 
Bernardi 2009), while galaxies in cluster centres may have 
had an unusual and extreme merger history (e.g. Boylan- 
Kolchin, Ma and Quataert 2006). Thus merger history and 
cluster environment may both influence the CMR. 

Our studies of the CMR in Roche, Bernardi and Hyde 
2009 (hereafter Paper I) suggested (see Section 3.1) that 
the radial colour gradients in E/S0 galaxies might be corre- 
lated with their other properties and it would be fruitful to 
investigate this further. These colour gradients are of inter- 
est in that they are sensitive tracers of evolution processes 
and vary strongly between the individual E/SOs. Most E/S0 
galaxies have negative colour gradients, meaning that the 
centres are reddest and colours become bluer outwards. This 
could be the result of a negative gradient in age, metallicity 
or dust, or some combination of these. 

Tamura et al. (2000) and Tamura and Ohta (2003) esti- 
mate the mean colour gradient in E/SOs as ^(tog"^) = — 0.09 
mag dex -1 (with a scatter of 0.04), and attributed this to 
a radial gradient of metallicity (rather than age) ^) = 
— 0.3±0.1, approximately constant to z ~ 1. Wu et al. (2005) 
similarly estimated d }^ g ^ = —0.25 with a large scatter 
a — 0.19. Mehlert et al. (2003) found negative metallicity 
gradients in Coma-cluster ellipticals, but negligible gradients 
in age or a/Fe. Michard (2005) concluded colour gradients in 
50 nearby ellipticals were primarily due to metallicity gradi- 
ents and the effects of dust were usually small. La Barbera 
and Carvalho (2009) compare optical and near-IR colour 
gradients and again conclude they are produced by negative 
gradients in metallicity (they find the radial age gradients 
in E/SOs may be even be positive). 

Kobayashi (2004) performed detailed chemodynamical 
simulations of a set of over 100 model elliptical galax- 
ies, with star-formation and merging, and predicted that 
spheroidals which formed monolithically, or at least with 
only minor mergers, would have steeper metallicity gradi- 
ents ( d^'iog V) — — 0-3 to -0.5) than those formed by major 
mergers ( ^ p° g f — —0.2). Within the simulated galaxy set, 
merger history was the primary determinant (more so than 
mass or luminosity) of present-day metallicity /colour gradi- 
ent. This might account for the observed wide distribution in 
colour gradients if similar numbers of spheroidals had mono- 
lithic, minor-merger or major-merger histories, and might 
also lead to an environmental dependence. 

In this paper we focus especially on Brightest Cluster 
Galaxies (BCGs) , of which several thousand have been iden- 
tified in the SDSS. Previously, the BCGs in the SDSS have 
been found to follow a steeper relation of radius to luminos- 
ity (r e ff oc L) than the other E/SOs (r e // oc L 0,6 ); (Bernardi 
et al. 2007; Bernardi 2009). In this paper, we look for sys- 
tematic differences in their CMR and radial colour gradients 



with respect to other E/S0 galaxies, especially those in the 
same luminosity range (M r < —22.5 to M r ~ —25). Some 
models, for example, predict a much flatter CMR for the 
BCGs, as a result of formation from a relatively large num- 
ber of mergers of (already) old and red galaxies (de Lucia 
and Blaizot 2007). Ko and Im (2005) found indications that 
the colour gradients of ellipticals in dense cluster environ- 
ments tended to be less strong than in field environments, 
which again could be due to (more) mergers. 

In Section 2 of this paper we describe the data used and 
our selection of early-type galaxies and BCGs. In Sections 3 
and 4 we investigate and compare their CMR and colour-a 
relation. In Section 5 we use an estimator of colour gradient 
based on the ratio of g and r-band effective radii to compare 
the colour gradients of BCGs and other E/SOs and examine 
the dependences on luminosity, radius, and density. In Sec- 
tion 6 we look at the influence of stellar age, as estimated 
from the spectra. We conclude in Section 7 with summary 
and further discussion. 

SDSS magnitudes are given in the AB system where 
rriAB = —48.60 — 2.5 log F„ (in ergs cm _2 s _1 Hz _1 ); equiv- 
alently, uiab = is 3631 Jy. We assume throughout a spa- 
tially flat cosmology with Ho = 70 km s _1 Mpc _1 , Qm = 
0.27 and S7a = 0.73, giving the age of the Universe as 13.88 
Gyr. 



2 SDSS DATA AND GALAXY SAMPLE 
SELECTION 

We first describe the selection of E/S0 and BCG galaxies 
used in the colour-magnitude relation analysis. For the fol- 
lowing analysis of colour gradients, we exclude some objects 
to give somewhat smaller samples, as described in Section 
5.1. 

In Paper I we selected 70378 E and SO galaxies out of a 
total of 367471 galaxies in the DR4 spectroscopic sample of 
the Sloan Digital Sky survey, with parameters updated to 
DR6. Strict criteria were applied; the E/SOs had to have de 
Vaucouleurs profile fractions (fracdev) from the SDSS disk- 
plus-bulge fits of 1.0 in both g and r (i.e. galaxies with any 
significant disk component are excluded), 'eclass' spectro- 
scopic classification parameters < (signifying absorption- 
line spectra), and dereddened r-band model magnitudes 
14.5 < r < 17.5. Redshift range is < z < 0.36. The 
same sample is used here, and (as in Paper I) we include the 
small 'sky subtraction' corrections from Hyde and Bernardi 
(2009); applied to the SDSS photometry these have the ef- 
fect of slightly increasing the effective radii and brightening 
the model magnitudes (in all passbands) for the galaxies of 
larger angular size. 

Here we make use of two samples of BCGs, both selected 
from the SDSS. The first are the 'C4 BCGs', a subsample 
of the first-ranked galaxies within the catalog of clusters de- 
tected using a 'C4' cluster-detection algorithm on the SDSS- 
DR2 data (Miller et al. 2005). Bernardi et al. (2007) selected 
a subset of 286 (of the original list of 748) BCGs as those 
with M r < —21, velocity dispersion (a) measurements and 
at least 10 other M r < —21 galaxies within 1.4 Mpc. We 
further restrict the Bernardi et al. (2007) C4 sample to the 
166 BCGs which lacked any spiral or disk component and 
had no neighbours close enough to confuse the radial profile. 
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These cover the redshift range z — 0.03142 to 0.16315 with 
a mean z = 0.0848. 

For the C4 BCGs (only) we have improved photometry 
available; the galaxies were individually re-examined using 
software developed with 'Matlab', which fitted new 'Gal- 
morph' magnitudes and radii in a system equivalent to the 
SDSS model magnitudes (Hyde and Bernardi 2009), with 
more accurate sky subtraction (not requiring further correc- 
tions) . 

The second, larger sample of BCGs is selected from the 
'max-BCG' catalog of Koester et al. (2007), and previously 
studied by Bernardi (2009); as in Bernardi (2009) the selec- 
tion criteria are to use only galaxies with SDSS spectra and 
which have frac<jcv > 0.8 in both g and r, 60 < a < 450 
km s~\ a concentration index > 2.4, and a dereddened 
model magnitude 14.5 < r < 17.5. This provides a sample of 
4919 BCGs covering a deeper redshift range z = 0.04144 to 
0.31722 with a mean z = 0.1837. For the max-BCG galaxies 
we have only the standard DR6 photometry and effective 
radii, to which we applied the sky-subtraction corrections 
from Hyde and Bernardi (2009). 

Fourthly, we select a comparison sample of high lumi- 
nosity non-BCG E/SOs. This is the subset of the 70378 E/S0 
galaxies with absolute magnitudes M r < —22.5 (there are 
20110) after excluding all galaxies identified as BCGs i.e. 
members of either the C4 BCG or max-BCG lists, to give a 
non-BCG sample of 18225. 



3 COLOUR-MAGNITUDE RELATION 
3.1 Background 

Bernardi et al. (2003 and 2005) investigated the CMR of 
early-type galaxies in the SDSS, deriving the rest-frame 
g — r colours from the observer-frame photometry using 
k-corrections calculated from models and template spectra 
(e.g. Hogg 2002). In Roche, Bernardi and Hyde (2009, here- 
after Paper I), we investigated the CMR of a larger sam- 
ple of 70378 E/S0 galaxies from later SDSS releases (Adel- 
man and McCarthy et al. 2006, 2008), again in terms of 
the mean rest-frame g — r colour as a function of red ab- 
solute magnitude (M r ). This time, rather than using mod- 
elled k-corrections, we obtained rest-frame colours using the 
observed SDSS spectra of the individual galaxies (indepen- 
dently of any models) , using two methods described briefly 
below (see Paper I for more details). 

Firstly, we extracted rest-frame g and r magnitudes 
directly from the spectra by integrating directly over the 
relevant wavelength ranges. We term these spectra-derived 
colours. This method does not use the imaging data at all 
or require any k-correction. However, the spectra-derived 
colour only measures g — r in the central 3 arcsec diameter 
aperture covered by the spectrograph fibre. As the angular 
size of an SDSS galaxy may be much larger than this, and 
will be correlated with luminosity, this can introduce a bias 
into the CMR (e.g. changing its slope) if the galaxies have 
radial colour gradients. 

The second method was to integrate the spectra over the 
observed and rest-frame passbands; the difference of these 
two magnitudes provided a spectra-derived k-correction for 
each individual galaxy. The SDSS imaging photometry pro- 
vides magnitudes based on profile fits to the galaxies (in 



the case of these galaxies, de Vaucouleurs profiles) , with the 
effective radius fixed across all passbands to that r e // mea- 
sured by fitting in the r-band. These are termed 'model' 
magnitudes and colours. Subtracting the spectra-derived k- 
corrections from the model magnitudes will then give rest- 
frame colours for the integrated emission from the whole of 
each galaxy. However, the model magnitude colours can also 
be subject to biases - in the opposite direction - because a 
k-correction based on an fixed-aperture spectrum is being 
applied to a colour derived from a model fit, which may be 
from a much larger aperture. 

In agreement with previous studies, we found approxi- 
mately linear CMRs, redder g — r colours for the more lumi- 
nous galaxies, and a mild blueward shift of the CMR with 
increasing redshift. However, our CMR based on spectra- 
derived colours was steeper than the CMR from model- 
magnitude colours ( d ^ M r ^ — —0.026 compared to -0.018), 
and showed more evolution with redshift. Similar variation 
in the CMR slopes between aperture and 'total' photomet- 
ric systems had been reported on by Scordeggio (2001) and 
Bernardi et al. (2003). The colour differences (g — r) spec — 
(g - r)model and (g - r) fib - (g - r) mode i will be sensitive to 
radial colour gradients. 

In Paper I we also examined the relation of colour to in- 
ternal stellar velocity dispersion a, the CaR. As expected, 
galaxies with greater a were redder, but we did not see a 
difference in slope between the spectra-derived and model- 
magnitude based CaR; both gave ^ ~ 0.19. This sug- 
gested that colour gradients might depend differently on lu- 
minosity and on a, and might have a positive correlation 
only with luminosity. We found that the difference between 
the 3-arcsec aperture and model magnitude colours was cor- 
related with the residual of a relative to the mean a-M r 
relation, and similarly with the M r residual to the M r -a 
relation. These two correlations indicated a stronger (nega- 
tive) colour gradient was associated with a higher luminosity 
and/or a lower a relative to the mean (a\M r ). We investi- 
gate this further in Section 5.2 below, after again looking at 
CMRs. 

In this paper we again make use of rest-frame g — r 
colours and absolute magnitudes M r obtained from the 
SDSS model magnitudes (or the Galmorph magnitudes in 
the case of the C4 BCGs), corrected to rest-frame using 
the spectra-derived k-corrections calculated as described in 
Paper I (including the 'signal/noise' correction, which had 
the effect of reducing the g-band k-correction by a few hun- 
dredths of a magnitude for the fainter objects, on the basis 
of an apparent discrepancy between the spectral calibration 
and the imaging photometry). We also make use of the rest- 
frame colours derived directly from the spectra. 

3.2 The CMRs of BCGs compared with other 
early-type galaxies 

Figure 1 shows the CMR, in model-magnitude g — r colours 
k-corrected to rest-frame and averaged in A(M r ) = 0.5 in- 
tervals, for the full E/S0 sample of Paper I (divided into 7 
redshift intervals), the C4 BCGs and the max-BCGs. As in 
Paper I the absolute magnitudes are 'de-evolved' by adding 
+0.86z magnitudes to separate the colour evolution from 
the effect of luminosity evolution. 

The C4 and max-BCG galaxies have similar CMRs, 
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Figure 1. Colour-magnitude relation, with model- magnitude 
colours averaged in A(M r ) = 0.5 intervals, of Brightest Cluster 
Galaxies (BCGs), separately for the lower-redshift C4 and the 
higher-redshift max-BCG samples, with the full sample of E/SOs 
divided by redshift. 



with the lower redshift C4s being slightly (~ 0.01 mag) red- 
der. Both sets of BCGs are on average slightly redder than 
the E/SOs at a similar luminosity and redshift and their 
CMRs are also flatter in slope. 

As in Paper, I we fit the galaxy CMRs with the form 



{(9 — r) r f\M r , z) = do + a-i M r + a 2 z 



(1) 



with separate linear dependencies on M r (01, the slope) and 
z (d2, the evolution) and no cross-term. The first line of Ta- 
ble 1 gives the fitted parameters for the (model-magnitudes) 
CMR of our full E/SO sample; these are the same as given 
in Paper I. 

Figure 2 compares the CMRs of the max-BCG and the 
M r < —22.5 non-BCG galaxies, divided into two intervals of 
redshift, and the C4 BCGs (in a single redshift interval - all 
are at z < 0.18); this is shown for both model-magnitude and 
spectra-derived colours. The two BCG samples have quite 
similar CMRs. In model magnitudes, the non-BCG CMR is 
significantly bluer than that of either BCG sample (in each 
redshift interval), except perhaps at the highest luminosities, 
and appears steeper in slope. In spectra-derived magnitudes, 
all the CMRs are moved slightly redwards (this could be the 
effect of negative colour gradients), but there is much less 
colour difference in between the BCGs and non-BCGs 

Table 1 gives the parameters fitted to these CMRs (an 
evolution term is only fitted for the max-BCG sample of 
BCGs as the C4 sample is too small and covers a narrow red- 
shift range). These indicate that the CMR of (both sets of) 
BCGs, although not entirely flat, has a slope (ai) less than 
half that of the non-BCG E/SOs of similar luminosity. This 
difference in a\ is seen for both model and spectra-derived 
colours. In model magnitudes we see little or no colour evo- 
lution (02), but in spectra-derived magnitudes there is (as 
noted in Paper I) some blueward evolution with redshift. 
This colour evolution is seen here for both BCGs and non- 
BCGs, with marginal evidence of a slightly lower rate for 
the BCGs. 

The mean residual of the max-BCG g~r, relative to to 
the (g - r\M r , z) fitted to the non-BCG M r < -22.5 E/SOs, 
is 0.0090±0.0007 in model magnitudes and 0.0016±0.0005 in 
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Figure 2. Colour-magnitude relation, in (a) model magnitudes, 
(b) spectra-derived magnitudes, for the max-BCGs, the sample 
of 18225 M r < -22.5 non-BCG E/SOs (divided into two intervals 
of redshift) and the C4 BCGs. 



spectra-derived magnitude. For the smaller C4 sample these 
A(g-r) are 0.0208±0.0030 and 0.0023±0.0029 magnitudes. 
Thus the BCGs are significantly redder than other E/SO 
(matched in luminosity) in their mean model-fit colours but 
with much less or no difference in their spectra-derived (i.e. 
central 3 arcsec aperture) colours. 
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Sample ao ai (slope) 02 (evolution) 

All E/SO (model mag) 0.3794 ± 0.0093 -(0.0174 ± 0.0005) -(0.0743 ± 0.0075) 

C4 BCGs (model mag) 0.6604 ± 0.1084 -(0.0056 ± 0.0047) 

max-BCGs (model mag) 0.5877 ± 0.0403 -(0.0082 ± 0.0018) -(0.0040 ± 0.0172) 

non-BCG (model mag) 0.3132 ± 0.0374 -(0.0200 ± 0.0017) -(0.0402 ± 0.0133) 

C4 BCG (spectra mag) 0.7121 ± 0.1028 -(0.0049 ± 0.0045) 

max-BCGs (spectra mag) 0.5118 ± 0.0282 -(0.0142 ± 0.0013) -(0.1825 ± 0.0121) 

non-BCG (spectra mag) 0.2234 ± 0.0286 -(0.0270 ± 0.0013) -(0.2326 ± 0.0104) 



Table 1. Fits to colour-magnitude relations 



Sample ao ai (slope) (evolution) 



All E/SO (model mag) 


0.3285 


± 


0.0044 


0.1914 


± 


0.0021 


-(0.1289 


± 


0.0053) 


C4 BCGs (model mag) 


0.5610 


± 


0.0859 


0.0946 


± 


0.0359 








max-BCGs (model mag) 


0.6584 


± 


0.0260 


0.0474 


± 


0.0109 


0.02049 ± 0.01528 


non-BCG (model mag) 


0.2270 


± 


0.0137 


0.2288 


± 


0.0059 


-(0.0664 


± 


0.0114) 


C4 BCG (spectra mag) 


0.6658 




: 0.082 


0.0656 


± 


0.0343 








max-BCGs (spectra mag) 


0.6538 


± 


0.0182 


0.0733 


± 


0.0077 


-(0.1375 


± 


0.0100) 


non-BCG (spectra mag) 


0.4381 


± 


0.0108 


0.1661 


± 


0.0046 


-(0.1985 


± 


0.0090) 



Table 2. Fits to colour-cr relations 
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Figure 3. Colour-cr relation with model-magnitude colours, of 
Brightest Cluster Galaxies (BCGs), separately for the lower- 
redshift C4 and the higher-redshift max-BCG samples, with the 
full sample of E/SOs divided by redshift. 



4 COLOUR- VELOCITY DISPERSION 
RELATION 

We can make the same comparison for the relation of colour 
to internal stellar velocity dispersion (cr), hereafter the CcrR. 
As described in Paper I, we apply small corrections to the 
SDSS a measurements to convert from the spectrograph 
aperture to a physical scale of r e ff /8, assuming the a oc r~ 0,04 
of Jorgensen et al. (1995). 

Figure 3 shows the colour-cr relation (CerR), defined in 
model magnitude rest-frame g — r, for the full E/SO sample 
and the two BCG samples. As in the CMR, we see the BCGs 
tend to be slightly redder (at a given a) and have a flatter 
CcrR. We fit with 



((<?• 



Iff 



\a, z) = ao + ai log a + <i2 z 



(2) 



The first line of Table 2 gives the fitted parameters for the 
whole E/SO sample (as given in Paper I). Figure 4 com- 
pares the CcrR of the two BCG samples with the non- 
BCG AI r < —22.5 galaxies, computed for both model and 
spectra-derived magnitudes. Table 2 gives the fitted param- 
eters (again, evolution is not fitted for the small C4 sam- 



ple). Again the BCG 



slope is less than half that of 



the non-BCG E/SOs. As we saw for the CMR, the spectra- 
derived colours show more blueward evolution than the 
model-magnitude colours, with the BCGs having a slightly 
lower evolution rate. 

In the CcrR the mean colour residual of the max-BCGs 
relative to the non-BCG (g-r\a, z) is 0.0070±0.0007 mag in 
model-magnitude g — r and 0.0016 ± 0.0005 mag in spectra- 
derived colour, and for the C4 sample these residuals are 
0.0187 ±0.0029 and 0.0042 ±0.0028 mag; very similar to the 
colour offsets in the CMR. 



5 COLOUR GRADIENTS IN E/SOS AND BCGS 

5.1 Method of Colour Gradient Analysis 

As a simple and rapidly computable quantifier of radial 
colour gradient, we make use of the ratio of the de Vau- 
couleurs effective radii fitted in the g band and the r band. 
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Figure 4. Colour-cr relation of the C4 BCGs and max-BCGs, 
compared with the 18225 M r < -22.5 non-BCG E/SOs. Shown 
for (a) model and (b)spectra-derived magnitudes. 



This ratio will be unity for a galaxy with no colour gradient 
and > 1 for a negative colour gradient, i.e. with the cen- 
tre redder than the outer annuli. A ratio of de Vaucouleurs 
r e ff (in V and K) was previously used as a colour gradient 
measure by Ko and Im (2005). 

To be redshift- independent, we need to use i e ff(g) and 
r e //(r) as they would be measured in the the rest-frame. 
We estimate these by linear interpolation between the r e ff 
fitted in the four passbands griz, using in the interpolation 
the passband mean wavelengths 4680, 6180, 7500, 8870A. 

We adopt a colour gradient measure as the restframe- 
corrected 



— 1. This will be related to but not directly 

equivalent to ^ 1 ^~ r r j , a measure of colour gradient obtained 
by directly fitting to the colour profiles of each galaxy. The 
relation between the two measures will depend on the exact 
form of the intensity profile and its variation with wave- 
length, and this may vary between individual ellipticals. To 
estimate this relation for the SDSS sample and to verify 



whether 



1 is a valid measure of colour gradient 



across this range of galaxies, we obtain from the SDSS Data 
Release 7 the tabulated g and r profiles for a small sample 
(90 galaxies). 

For the 90 galaxies, which include 21 BCGs, we make a 
linear error-weighted fit of the form 



(g — r) — ao ± di (log radius) 



(3) 



Colour Gradients of BCGs and E/SO Galaxies 7 



0.2 



-0.2 



S -0-4 



-0.6 



-0.8 



All galaxie 

r e „>1.8 


s correlation 
arcsec only 


-0.0119-0.9398*x 

-0.0089-0.8736*x - 

r eff <1 .8 arcsec * 
r eff >1.8 arcsec • 


'■- x * 

\ % 

K * * 


.* ■ 
■ 


J f 















0.1 0.2 0.3 0.4 

Our colour gradient measure [r 8 ff(g)/r eff (r)]-1 



0.5 



Figure 5. The relation between two measures of colour gradi- 
ent for an early-type galaxy, the ratio of de Vaucouleurs radii 
Iei f , 9 , — 1, used in this paper, and 4rf — ~r from a fit to the 

r eff\ r ) "(tog r) 

colour profile, for 90 galaxies from our sample with symbol types 
indicating apparent size (in the r band). Best-fit linear correla- 
tions are plotted for all 90 and for the 51 galaxies with r e f t > 1.8 
arcsec. 
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Figure 6. As Figure 5 but with symbols indicating that galaxies 
are BCGs, higher luminosity (M r < -21.5) but non-BCG E/SOs, 
or lower lminosity. Best-fit linear correlations are shown for each 
of these three classes. These correlations are quite similar but the 
scatter is less for the BCGs. 



to the g — r colour profile over all annular bins from an inner 
radius 0.67 arcsec to the outer limit of the tabulated profile, 
usually 28 arcsec, examining each fit by eye. 

Figure 5 shows the gradient ^f ~ r ^ plotted against our 



colour gradient 



ettM 



1. It is clearly correlated with 



a linear correlation coefficient 0.768 and best fit relation 
-0.0119(±0.0038)-0.9398(±0.0336)a;, passing very close to 
the origin as expected, but the scatter of 0.122 is rather 
large. This is in part because a significant minority of the 
ellipticals have colour profiles not linear in log radius; we see 
non-monotonic gradients for about 10% and in these cases, 
the r e f f ratio tends to be sensitive to the colour profile at a 
larger radius than the direct fit (which is weighted towards 
smaller radii where the error bars are smaller). 

Despite this, BCGs and high and low luminosity E/SOs, 
considered separately (Figure 6), show similar relations be- 
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Figure 7. Residuals of fitted gradient 

lation -0.0119 - 0.9398[ ^ e// ^| - 1], plotted against (a) absolute 
magnitude M r and (b) apparent size r e ff. 



tween the two colour gradient measures. However, the re- 
lation is much tighter for the BCGs with a scatter of only 
0.053, compared to 0.123 (high L E/S0) and 0.150 (low L 
E/S0). As a further test, we examine (Figure 7) the residu- 



als of the 



to the best-fit linear relation, as a function 



liflog r) 

of M r and apparent r e ff (in the r-band); a significant trend 
could introduce a bias into the investigations of this paper. 

We find the scatter tends to increase at lower luminosity 
but there is not an obvious systematic trend. Neither do we 
see a trend with apparent radius, down to about 1.8 arcsec. 
However, for galaxies smaller than this radius the scatter 
increases greatly and asymmetrically, which is probably un- 
acceptable for our purposes. For galaxies with r e // < 1.8 
arcsec it may not be possible to meaningfully characterize 
the colour gradient at all using SDSS data. For early-type 
galaxies of larger apparent size, including the BCGs, the 
ratio of de Vaucoleurs radii does appear to be a valid, con- 
sistent and approximately linear colour gradient measure. 

Returning to Figure 5, if all r e // < 1.8 arcsec galaxies 
are excluded, leaving 51, the best-fit linear relation changes 
only slightly to -0.0089(±0.0039) - 0.8736(±0.0346)a;, but 
the scatter is decreased by more than a factor of two to a 
much more acceptable 0.056 (similar to the scatter for BCGs 
only), and the correlation coefficient is increased to 0.862. 

Some galaxies in our sample have very strong (> 0.5) 
or inverted (< 0) colour gradients (i.e. r e //(g) < r e //(r)). 
Examining some of these, we find that some are galaxies 
confused with or very close to (red or blue) stars, and that 
there are also some misclassified spirals present, despite our 
strict selection criterion of fracDev = 1 in both g and r. To 
reduce any biases that these types of objects may introduce, 
we exclude some galaxies from the colour gradient analysis. 
Firstly, we excluded the objects with gradients (i.e. 
"f f ^ — 1) calculated as < —0.5 or > 1, far out of the 



oft 



normal range (there were only 130/70378). Secondly, we ex- 
amined by eye all galaxies in the E/S0 sample with gradients 
0.4-1.0; these numbered 1869 out of 70378. It was estimated 
that 108 of these were elliptical-elliptical mergers, 26 were 
mergers involving spirals, 98 were barred spirals, 142 were 
edge-on disks, 120 were other galaxies with visible spiral 
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arms and 88 were confused with stars. All of these except 
the elliptical mergers, 474 galaxies, were marked for exclu- 
sion. It should also be noted that many other high colour 
gradient galaxies were asymmetric or disturbed. 

At lower colour gradients the proportion of face-on spi- 
rals decreases steeply (13% at gradient > 0.4, 5% at 0.3, 
2% at 0.2) but less steeply for edge-on disks (8% at gradient 
> 0.4, 8% at 0.3, 3% at 0.2). However, the latter could easily 
be removed as they had axis ratios (from the r-band de Vau- 
couleurs fits) of typically b/a ~ 0.3 and always b/a < 0.4, 
whereas virtually all ellipticals had b/a > 0.4. We therefore 
exclude all galaxies with b/a < 0.4, whatever their colour 
gradient. Lastly, because of their very large scatter on Fig- 
ure 7, we exclude all objects with r-band r e // smaller than 
1.8 arcsec. 

These exclusion criteria, for the purposes of colour gra- 
dient analysis, reduce the all E/S0 sample to 37001, the 
M — r < —22.5 non-BCG sample to 13652, and the max- 
BCG sample to 4485 (we note that these exclusions make 
very little difference to the colour-magnitude relations of 
the respective samples, the main effect being increased er- 
ror bars at the lowest luminosities). 



5.2 Results 

After this filtering, for the whole E/S0 sample the mean 



1 is 0.11115±0.00061 with a scatter 0.11724. From 



the relation on Figure 6 this corresponds to d( - 9 r - > 



d(log r) 

and so lies half-way between the mean gradient t ^ 3 ~ r - ) 



0.10 



<2(log r) 

0.071 measured for similarly selected E/SOs by La Barbera 



-0.14 found 



and Carvalho (2009) and the larger d( - log ^ 
by Suh et al. (2010) for 'red-cored' morphological early-types 
after the exclusion of zero and inverted gradients. 

For the M r < —22.5 non-BCGs sample the mean gradi- 
ent is similar, 0.10550 ± 0.00098 with a scatter 0.11414. But 
for the max-BCGs the mean gradient is significantly lower 
at 0.08142 ± 0.00114 with scatter 0.07645. Hence, we find 
evidence that BCGs have flatter colour gradients than other 
E/SOs; the difference is modest but with this large sample 
highly significant - the BCGs' mean gradient is 22.8 ± 1.8% 
less than for the M r < -22.5 non-BCGs and 26.7 ± 1.8% 
less than for all E/S0. 



Figure 8 shows the mean 



1 as a function of 



absolute magnitude M r . The luminosity dependence is mild. 
There is a broad maximum at M r ~ —22, with a decrease 
of up to 20% to higher luminosity (BCGs excluded) and 
24% to lower luminosities. BCGs have a significantly lower 
mean colour gradient than other E/SOs across their range of 



luminosities. Figure 9 shows the mean 



e££M 



1 with the 



galaxies divided by velocity dispersion a. Colour gradient is 
maximum at low velocity dispersions of 100-150 km s _1 and 
decreases by almost 1/2 to a ~ 300 km s _1 . The BCGs have 
a lower colour gradient that non-BCGs in each a interval, 
although the gap between the two is narrower (~ 0.01) than 
with the galaxies in M r intervals (0.02-0.03). In the non- 
BCG sample, which has a sharp lower cutoff in luminosity, 
there is a steep increase in colour gradient at a < 180 km 
s _1 . This reflects the finding in Paper I that a low a relative 
to luminosity is associated with a strong colour gradient. 
This is again seen on Figure 10 where the E/S0 sample 
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Figure 8. The mean 



e//(g) 



1 (corrected to rest-frame), a 



measure of colour gradient, for the full E/S0 sample, the M r < 
—22.5 non-BCGs and the max-BCGs, as a function of absolute 
magnitude M r . 
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a measure of colour gradient, for the full E/S0 sample, the M r < 
—22.5 non-BCGs and the max-BCGs, as a function of velocity 
dispersion a. 



is divided by absolute magnitude. Strong (> 0.15) colour 
gradients are a feature of galaxies with moderate or high 
luminosity of M r < —21 combined with log a < 2.25 [p < 
180 km s _1 ). 

Figure 11 shows colour gradient against effective radius. 
For the E/SOs the mean ( te ff, . — 1 almost doubles from the 
smallest radii up to r e ff = 8-10 kpc, and similarly for the 
M r < —22.5 non-BCGs, for which the gradient peaks at 



larger radii of 11-14 kpc (presumably because of the greater 
mean luminosity of this sample). In both samples, mean 
gradient decreases at even larger radii. In BCGs the colour 
gradient shows a much weaker variation with radius and 
remains lower than in non-BCGs, especially at 5 < r e g < 20 
kpc. 

To better understand the non-monotonic variation of 



1 with r c ff (r), and remove possible selection effects 
(from the flux limit of the sample) we re-examine this rela- 
tion with the E/SOs divided by M r (Figure 12). For mod- 
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Figure 12. The mean Jp: — 1 (corrected to rest-frame) for 

the full E/SO sample divided by M r , and for the max-BCGs, as 
a function of red-band effective radius r e ff . 
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Figure 11. The mean ^ e ^|^j — 1 (corrected to rest-frame), a 
measure of colour gradient, as a function of red-band effective 
radius r e ff for the full E/SO sample, the M r < —22.5 non-BCGs 
and the max-BCGs. 



erate luminosities there is simply a steep increase in colour 
gradient with radius. The turnover at large radii is produced 
by M r < —23 galaxies, in which colour gradient is 'sup- 
pressed' and the correlation with radius is flattened. This 
effect is also seen, even more strongly, in the BCGs. 

In Paper I, by comparing model and aperture colours, 
we found evidence that stronger colour gradients were asso- 
ciated with E/SO galaxies with high luminosity relative to 
a. Figure 13 plots ^ e ^|^ — 1 against 10 log a + M r (approx- 
imately the a residual relative to the a-M r relation; Forbes 
and Ponman 1999). In the full E/SO sample the colour gra- 
dient increases by a factor of two from the highest to the 
lowest ratios of a to luminosity. In the M r < —22.5 non- 
BCG sample this increase is an even larger factor of three. 
In marked contrast, the colour gradient in BCGs shows little 
or no dependence on the a residual and remains at ~ 0.08. 

Figure 14 shows color gradient against a 2 /r^ s , a more 
physical quantity as it is a direct measure of mass density 
(subject to a mild dependence on long-axis orientation). We 
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Figure 13. The mean **** ^ — 1 (corrected to rest-frame), for 
the full E/SO sample, the M r < -22.5 non-BCGs, and the max- 
BCGs, as a function of 10 log a + M r . 

express this in the form log (c 2 /rg ff ) for a in km s _1 and 
r e ff in kpc. In these units the mean density for the E/SOs 
is 3.38 with a scatter 0.42. This can also be converted to 
units of solar masses Mq. From Boylan-Kolchin, Ma and 
Quataert (2006) the mass within r < r e ff is approximately 
2.90 " g eff , where the gravitational constant G = 4.30 x 10~ 3 
pc Mq km s -1 ) 2 . Within r e ff the volume is Im^g, and the 
mean mass density, log (M Q /kpc 3 ) = log (cr 2 /r 2 ff ) + 5.21. 
The full E/SO sample and the M r < -22.5 non-BCGs show 
a similar decrease in colour gradient with increasing mass 
density, whereas for BCGs r " ff( f\ — 1 ~ 0.08 over the whole 

r e ff(r) 

density range. 

Figure 15 shows colour gradient against density for 
E/SOs divided by absolute magnitude. The negative correla- 
tion with density is seen over a very wide luminosity range, 
and strongest at —21 > M r > — 23. 

To briefly summarize at this point, we find colour gra- 
dients in E/SOs to (i) be greatest at intermediate luminosi- 
ties and (ii) at a fixed luminosity to be anticorrelated with 
velocity dispersion and mass density, and that (iii) colour 



10 N. Roche, M. Bernardi and J. Hyde 



All E/SO 
non-BCG M r <-22.5 
max-BCGs 




1.5 



2.5 



3.5 



4 



4.5 



log Density (sigma[km/s] / r[kpc]) 



Figure 14. The mean 



ef/(g) 



1 (corrected to rest-frame), for 
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pressed in a log scale for units of km s~ 1 and kpc. 
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by M r , and for the max-BCGs, as a function of a / r^ xp , 
of mass density. 



a measure 



gradients tend to be flatter in BCGs and to show much less 
dependence on galaxy properties. We now investigate the 
effects of age. 



6 AGE EFFECTS ON COLOUR GRADIENT 
AND COLOUR 

For most galaxies in the E/SO sample (64774/70378) and 
the M r < -22.5 non-BCG sample (16337/18225), we have 
mean stellar age ('age50') estimated from the spectral in- 
dex analysis of Gallazzi et al. (2005) performed on the 
3.7 x 10 5 galaxies of the DR4. This catalog includes only 
part of the max-BCG list. By position-matching the two 
lists we find there are stellar age estimates for almost half 
(2274/4919) the max-BCG sample, sufficient to investigate 
age-dependent effects. These stellar ages will generally be 
less than the lookback time to galaxy formation, and as they 
are luminosity-weighted, the age estimate of a galaxy may 



be greatly reduced if it has experienced recent episodes of 
star-formation. 

The stellar mean age (and scatter) are 6.92 (1.74) Gyr 
for all E/SOs, 7.04 (1.69) Gyr for the non-BCGs, 7.54 (1.39) 
Gyr for the max-BCGs and 8.12 (1.33) Gyr for the C4 BCGs. 
The greater mean age of the C4 BCGs is due to these being 
a lower redshift sample than the others. Figure 16a shows 
the age distributions for the 3 deep samples; the modal ages 
are similar but there are fewer younger (< 6 Gyr) and more 
older ages (9-10 Gyr) ages amongst the BCGs. 

We add the mean stellar ages to the lookback times at 
the redshifts of observation, giving lookback times to the 
mean stellar formation time, and then convert these back to 
'formation' redshifts (not the redshift of galaxy formation 
but the mean redshift of luminosity- weighted star formation, 
which will be lower) for the galaxies. Figure 16b shows the 
distribution of the mean redshift of star formation; there is 
more of a difference here between the all-E/SO sample and 
the BCGs and other high luminosity galaxies, as the more 
luminous galaxies tend to be observed at higher redshifts 
as well as being older. The mean star-formation redshift is 
1.32 for all E/SO, 1.56 for the non-BCGs, 1.53 Gyr for the 
C4 BCGs and 1.82 for the max-BCGs; mean lookback times 
are 8.47, 9.17, 9.22 and 9.85 Gyr. 

Figure 17 shows mean age against M r \ the BCGs are 
not only older than other E/SOs but show less of an increase 
in age with luminosity. 

We investigate whether the redder colours and weaker 
colour gradients of BCGs are related to their greater ages, 
and in addition, whether the trends in (non-BCG) colour 
gradients with galaxy properties are caused by age effects. 
A higher luminosity or lower a relative to the mean (a\M r ) 
relation, found (in Paper I and here) to be associated with 
a stronger colour gradient, is also correlated with a younger 
stellar age (Forbes and Ponman 1999; Bernardi et al. 2005; 
Gallazzi et al. 2006). 

On Figure 18, for non-BCG E/SOs the colour gradients 
are maximum at stellar ages 4-5 Gyr and decrease by almost 
a factor of two to the highest ages. On Figure 19 it can be 
seen that this peak in gradient vs. age is most pronounced 
for galaxies of intermediate luminosities of —21 > M r > 
-23, weaker at both -20 > M r > -21 and -23 > M r > 
—24, and apparently absent at M r > — 20 and M r < —24, 
following a rather symmetric luminosity trend. The BCG 
colour gradients show no age dependence (at > 4 Gyr), and 
again the BCGs have lower colour gradients than the non- 
BCG E/SO sample, at all but the highest ages where the two 
converge. 

A potential problem with this analysis is that the stellar 
ages are 'central' ages estimated from the spectra in fixed 
apertures, and if colour gradients were typically driven by 
age gradients this would bias the age-gradient relation. We 
are not able to investigate age gradients with the single aper- 
ture of SDSS spectroscopy, but Tamura and Ohta (2004) and 
La Barbera and Carvalho (2009) find that age gradients in 
early-types are typically very small or consistent with zero 
and the negative colour gradient is essentially a metallicity 
gradient (positive colour gradients can be stellar age related, 
as discussed below, but these occur in < 8% of these sam- 
ples) . 

At ages of < 3 Gyr there is sharp drop in mean colour 
gradient in the youngward direction, seemingly for all lu- 
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Figure 16. Distribution of the (a) mean stellar age (b) mean 
redshift of star formation, for the full E/SO sample, the M r < 
—22.5 non-BCGs, and the max-BCGs, normalized to the same 
area under the curves. 

minosities of E/SO and even the BCGs (although very few 
- 1.1% - of this BCG sample are this young). This can be 
explained by the younger stellar populations within these 
galaxies being more centrally concentrated, giving positive 
age gradients which counteract negative metallicity gradi- 
ents to flatten the colour gradients. Some post-starburst 
galaxies have significantly opposite gradients i.e. blue cores 
(Yamauchi and Goto 2005) from central populations of very 
young stars (Suh et al. 2010). 

We examined the images and spectra of some of the 
galaxies with very inverted colour gradients. Although some 
of the most extreme examples were galaxies confused with 
stellar images, out of 40 galaxies with gradient < —0.1 
and age < 3 Gyr, most looked like normal ellipticals and 
many had strong Balmer absorption lines indicative of 
post-starburst galaxies. The best 7 examples had equivalent 
widths of 5-8A for both H/3 and H<5 absorption (these are 
J034453.17-065442.9, J150020.11+040000.2, J211348.16- 
063059.6, J210258.87+103300.6, J225226.07+150303.1, 
J130029.44+545503.9 and J155435.53+291319.9); there 
were a further 6 with H/3 EW = 4.5-5.0 A, and two more 
with less absorption in H/3 but Hr5 EW > 5A. 

Figures 18 and 19 show that differences in the colour 
gradient properties of BCGs and non-BCGs are still appar- 
ent when they are age- matched, and so are not simply caused 
by the difference in their mean age. In the same way we can 
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Figure 17. The mean stellar age against absolute magnitude 
M r for the full E/SO sample, the M r < -22.5 non-BCGs, and 
the max-BCGs, and the less deep C4 BCG sample. 
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Figure 18. The mean colour gradient r"j^(r) — ^ (corrected to 
rest-frame) for the full E/SO sample, the M r < -22.5 non-BCGs, 
and the max-BCGs, as a function of mean stellar age ('agc50') 
in Gyr, for all galaxies for which we have age estimates from the 
methods of Gallazzi et al. (2005, 2006). 

investigate whether this is the case for their colours. Figure 
20a shows mean model g — r against stellar age; the youngest 
galaxies are much bluer, but the age dependence flattens at 
> 6 Gyr, while the BCGs are redder than age-matched non- 
BCGs. The M r < -22.5 non-BCGs are also redder than 
the full E/SO sample, due to the higher metallicity of the 
more luminous galaxies (hence these are the more appro- 
priate sample to compare colours with the BCGs), but the 
BCGs are ~ 0.01 mag redder still. Figure 20b shows how 
colour depends strongly on luminosity even when ages are 
matched (hence the CMR slope is not just age driven). 

The spectra-derived colours (Figure 20c) follow a more 
linear relation to stellar age with A(g — r) ~ 0.5A(age), 
and, notably, do not appear to differ between BCGs and 
age-matched, similarly luminous non-BCGs (except for the 
few very young BCGs, which are redder). The difference 
between the model and spectra-derived (aperture) colours 
is caused by colour gradients. 

We found BCGs to be at least A(g — r) ~ 0.01 mag 
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Figure 19. The mean colour gradient 



4 5 6 7 8 
Mean stellar age (Gyr) 

//(g) 



eZ/M 



10 11 12 



1 (corrected to 



rest-frame) for the full E/S0 sample, as a function of mean stellar 
age ('age50') in Gyr, with galaxies divided by absolute magnitude 

M r . 



redder than other E/SOs in the model-magnitude CMR and 
CaR (Figures 1, 2, 3 and 4) with a smaller (< 0.005 mag) 
or no significant difference in the spectra-derived colours. 
This can now be explained as a combination of (i) the 0.5 
Gyr greater mean age of BCGs compared with luminosity- 
matched non-BCGs, corresponding to A(g — r) ~ 0.0025 in 
the aperture colour, and (ii) the flatter colour gradients of 
BCGs, meaning that even if their spectra-derived (central) 
colours match those of non-BCGs, their model-magnitude 
colours will be redder. This could account for the remaining 
~ 0.0075 magnitudes. 

To further investigate the relationship between age and 

e//(g) 



colour gradient, we show mean 



e//M 



1 against radius 



(Figure 21) 10 log a + M r (Figure 22) and density (Figure 
23) for the full E/S0 sample divided by mean stellar 'for- 
mation' redshift. BCGs are included to see if their different 
properties with respect to other E/SOs resemble the effects 
of increased age. We also (Figure 24) plot colour gradient 

2 

against dynamic mass, estimated as 5.8^^^— ; for units of 
kms -1 and kpc, log (Md yn /M e ) = log r e ff + 21ogcr + 6.13. 

Colour gradients show a broad peak at Md yn — 
10 114 M Q , which is more pronounced for galaxies with 
younger stellar ages. Colour gradients in the galaxies with 
the youngest stellar populations, Zf orm < 0.5, show correla- 
tions with galaxy properties, similar to the 0.5 < zf orm < 
1.0 galaxies, but offset downwards, as might result from the 
presence of an admixture of post-starbursts with negative 
gradients. Colour gradients in E/SOs depend significantly on 
radius, density, mass etc. within each zf orm interval. How- 
ever, these correlations, especially those with radius and 
density, become less strong in Zf orrn > 1.5 galaxies. We dis- 
cuss this further in the next Section. 




4 5 6 7 
Mean stellar age (Gyr) 



0.84 
0.82 

0.8 
0.78 
0.76 - 
0.74 
0.72 - 

0.7 
0.68 
0.66 
0.64 
0.62 

0.6 

0.86 
0.84 
0.82 

0.8 
0.78 
0.76 
0.74 
0.72 - 

0.7 
0.68 



(b) 




E/S0 M r >-20 
E/S0 -20>M r >-21 
E/S0 -21>M r >-22 
E/S0 -22>M r >-23 
E/S0 -23>M r >-24 

E/S0 M r <-24 
max-BCGs 



4 5 6 7 8 
Mean stellar age (Gyr) 



10 11 12 



(c) 




All E/S0 
non-BCG M r <-22.5 
max-BCGs 





0.88 
0.86 
0.84 
0.82 

0.8 
0.78 
0.76 
0.74 
0.72 - 

0.7 

0.68 

0.66 

0.64 

0.62 




4 5 6 7 8 
Mean stellar age (Gyr) 



10 11 12 



(d) 




E/S0 M r >-20 
E/S0 -20>M r >-21 
E/S0 -21>M r >-22 
E/S0 -22>M r >-23 
E/S0 -23>M r >-24 

E/S0 M r <-24 
max-BCGs 



4 5 6 7 8 
Mean stellar age (Gyr) 



10 11 12 



Figure 20. The mean rest-frame g—r from the k-corrected model 
magnitudes as a function of mean stellar age (age50) from Gallazzi 
et al. (2005), (a) for the full E/S0 sample, the M r < -22.5 non- 
BCGs, and the max-BCGs, (b) for E/SOs divided by luminosity; 
(c,d) as (a,b) using spectra-derived magnitudes. 
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Figure 21. The mean colour gradient ^ — 1 (corrected to 

rest-frame) for the full E/SO sample divided by mean stellar age 
(converted to a mean stellar formation redshift), and for the max- 
BCGs, as a function of radius. 
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Figure 22. As Figure 21 as a function of the velocity dispersion 
to luminosity ratio. 
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Figure 23. As Figure 21 as a function of mass density. 
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(5.8r eff a 2 /G). 

7 SUMMARY AND DISCUSSION 

(i) The colour-magnitude relation (CMR) and colour- 
velocity dispersion relation (CcrR) of BCGs are significantly 
flatter in slope( aim jiogg ) than the respective rela- 

tions for non-BCG E/SO galaxies of comparable (high) lu- 
minosity. The difference in slope is about a factor of two, 
whether we define the g — r colour using k-corrected model 
magnitudes or rest-frame magnitudes derived from integrat- 
ing the SDSS spectra. BCGs are also, on average, ~ 0.01 
mag redder in model magnitudes (k-corrected to rest-frame) 
than E/SOs of the same M r or a, but with less or no offset 
in spectra-derived colours. 

The hierarchical merging model of Lucia and Blaizot 
(2007) predicts a flat CMR for BCGs as a result of late as- 
sembly from a large number of red progenitors which formed 
their stars very much earlier (z ~ 5). This model included 
a very early truncation of star-formation. Skelton, Bell and 
Somerville (2009) present a simplified model in which dry 
mergers of already merged E/SOs on a 'creation red se- 
quence' mildly flatten the CMR slope at higher luminosi- 
ties. There may indeed be some flattening in the bright end 
of the CMR for all E/SOs (see Paper I), but for BCGs we 
find much more, although the BCG CMR is not entirely flat. 
This suggests the BCGs or their progenitors evolve by a less 
extreme scenario than the first of these models but perhaps 
with a higher dry merger rate compared to the E/SOs in the 
second model. In addition to being flattened, the BCG CMR 
is offset redwards. 

(ii) As a simple quantifier of radial colour gradient we 
use 1°^^ ~ 1) with the effective radii r e //(p) and r e //(r), 
corrected to the rest-frame by linearly interpolating between 
the radii fitted in the bands g, r, i and z. We show that at 
least for galaxies with angular r e // radii above 1.8 arcsec, 
this ratio is well and linearly correlated to ^n^rjj across the 
observed range of colour gradients and galaxy luminosities, 
and estimate d ,a~ r \ ^ -0.87( w/( ,^ - 1). The large scatter 

"(log r ) r c//( r ) ' 

in this correlation, which we initially estimated as 0.122, is 
reduced to a much more acceptable 0.056 if only galaxies 
with r e ff larger than 1.8 arcsec are included. 

After excluding the r e // < 1.8 arcsec galaxies and a 
few others which we believed to be spiral types, we found 
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for the E/SOs a mean colour gradient, using our measure, of 
0.11115 ± 0.00061 with a comparable galaxy-to galaxy scat- 
ter 0.1172. For the M r < -22.5 non-BCG the mean is simi- 
lar at 0.10550 ±0.00098 (scatter 0.1141) and for max-BCGs 
the mean is 0.08142 ± 0.00114 (scatter 0.0764). Therefore 
we find at high significance that BCGs have a flatter (by 
an average of 23%) radial colour gradient than other high 
luminosity spheroidals. 

This again could be due to intensified merging in the 
cluster-core environment. La Barbera et al. (2005), and sim- 
ilarly Ko and Im (2005), reported mean colour gradients in 
E/SOs of all luminosities were weaker by almost a factor two 
in the richest clusters (N ga i > 50) compared to the field, 
and suggested this was because formation histories in denser 
cluster environments included more elliptical-elliptical merg- 
ers (e.g. Tran et al. 2005), which reduce colour/metallicity 
gradients (Kobayashi 2004). 

The simulations of Boylan-Kolchin, Ma and Quataert 
(2006) account for the steeper R e -L of BCGs, i.e. their 'ab- 
normally' large sizes (e.g. Bernardi 2009), as the result of re- 
peated 'dry' (dissipationless, non-star-forming) mergers, and 
especially of the anisotropic accretion of smaller spheroidals 
in near-radial, low angular momentum orbits, which might 
occur preferentially in a cluster core. A few BCGs have been 
observed undergoing such mergers (Liu et al. 2008; Tran 
et al. 2008). Di Matteo et al. (2009) predict from simula- 
tions that the metallicity gradient in a dissipationless ('dry') 
merger remnant will consistently be ~ 0.6 the mean of the 
progenitor galaxies' gradients (but not reduced to zero). 

(iii) In non-BCG E/SOs, we examine the trends in mean 
colour gradient as a function of other galaxy properties. We 
find a dependence of colour gradient on absolute magnitude 
M r , with a broad peak in colour gradient at M r ~ —22 and 
a decrease to lower and higher luminosities. The luminosity 
dependence is relatively mild, less than a factor 1.5, which 
may explain why it was not seen by La Barbera et al. (2005). 
However, Spolaor et al. (2009) do find a luminosity depen- 
dence similar to ours. Colour gradients tend to decrease with 
increasing velocity dispersion <r, by almost 1/2 between 150 
and 300 km s~ , confirming the suggestion in our Paper I 
that the trends with M r and a are different. 

(iv) Colour gradients in E/SOs show strong correla- 
tions with some other galaxy properties. Colour gradients 
increase, by about a factor of two, with from small to larger 
effective radius, up to a maximum at 8-12 kpc (depending on 
luminosity). A positive correlation with radius was prevously 
reported by Tamura and Ohta (2003) for a small sample of 
E/SOs in Abell 2199. At even larger radii we find the mean 
colour gradients decrease, but this is actually because of a 
reduced gradient in the most luminous (M r < —23) galax- 
ies. If E/SOs are divided by luminosity, within low/moderate 
luminosity intervals there is simply a steep increase in colour 
gradient with r c ff . 

We find colour gradients to be negatively correlated 
with 10 log a + M r (the a residual relative to (a\M r )) and 
mass density (a 2 /r 2 cB ). Of course, these quantities are re- 
lated - a high density implies a high a at a given dynamic 
mass (rather than luminosity). These negative correlations 
do not seem to be caused by a selection effect (from the flux 
limit of the sample), as they are seen within a wide range 



of narrow luminosity intervals, being strongest for E/SOs of 
moderate/high luminosity —21 < M r < —24. 

(v) We also examine the trend in colour gradient with 
the 'mean stellar ages' estimated (Gallazzi et al. 2005, 
2006) for most of these galaxies using 5 indices from the 
SDSS spectra. In non-BCG E/SOs the colour gradients are 
strongest for relatively young ages of 4-5 Gyr (which may 
signify the presence of younger stellar populations or more 
extended star- formation, rather than the galaxies being 
younger) and decrease by up to a factor of two to the old- 
est ages of 11 Gyr. This may be primarily because older 
stellar populations have experienced more mergers and in- 
teractions, especially at high redshifts. Early-type galaxies 
which ceased forming stars and joined the red sequence more 
recently (e.g. post-mergers of spirals at z < 1) would tend to 
have stronger colour gradients, if it is dry mergers (occurring 
after the end of star-formation) that flatten the gradients. 

However, there is a drop in mean colour gradient for the 
galaxies with mean stellar ages less than 3-4 Gyr (even in 
BCGs). This is likely due to the inclusion of post-starburst 
galaxies, which can have blue cores giving inverted (positive) 
colour gradients (Friaca and Terlevich 1999; Yamauchi and 
Goto 2005) together with young age indicators like strong 
Balmer absorption lines. To confirm this we examined the 
spectra of 40 galaxies with age < 3 Gyr and inverted colour 
gradients (< —0.1) and found many with enhanced Balmer 
absorption including 7 with equivalent widths > 5A in both 
H/3 and H5. Suh et al. (2010) found that blue-cored early- 
types in general had at least mildly enhanced H/3 absorp- 
tion, and showed that the Balmer enhancement and the 
positive colour gradients could both be accounted for by 
centrally concentrated, very young stellar populations (age 
< 0.5Gyr), comprising only ~ 0.5-2% of the stellar mass. 

(vi) Colour gradients in BCGs are consistently lower 
than in non-BCG spheroidals of the same luminosity, veloc- 
ity dispersion or radius. Some classes of non-BCG ellipticals, 
those with (i) the highest luminosities and largest radii, (ii) 
the highest a relative to M r , or (iii) the highest densities, 
also have mean colour gradients of r "^ g ~ 0.08. But within 

the BCG class mean T 'f^ 3 \ _ 1 remains at ~ 0.08, almost 
regardless of the other galaxy properties. The mean colour 
gradient in BCGs decreases a little for the largest radii and 
with increasing dynamic mass, but is uncorrelated with mass 
density, the 10 log a + M r residual, or stellar age (except for 
the drop at < 3 Gyr). 

The process of colour gradient flattening through dry 
mergers appears to have operated most strongly on (all) 
BCGs, reducing the mean colour gradients to ~ 0.08, but 
never much lower than this, whatever the initial value. In 
the Di Matteo et al. (2009) model of dry mergers, a galaxy 
with a strong colour gradient will probably merge with one 
of lower gradient and experience a significant flattening, but 
a galaxy with already a very flat colour gradient would prob- 
ably merge with one with higher gradient and experience no 
further decrease or even a slight increase. With repeated 
mergers this would cause an exponential decrease in the 
mean gradient and also a narrowed distribution, which is 
what we seem to find for the merged classes of galaxies, 
with mean gradients of ~ 0.08 consistently. We found the 
BCGs had | less galaxy-to-galaxy variation in colour gra- 
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dient than the E/SOs (this difference might be even greater 
after taking measurement errors into account). 

In addition, the radial mergers thought to occur in 
BCGs would strongly increase their radii (and reduce den- 
sity) while simultanously flattening colour gradients, thus 
erasing the gradient-radius correlation (and gradient-density 
anticorrelation) seen in other spheroidals. 

(vii) We examine colour as a function of stellar age, 
comparing the BCGs and other E/SOs. The BCGs are on 
average 0.5 Gyr older at observation than the M r < —22.5 
non-BCGs, and their mean lookback time to formation is 
estimated as 0.68 Gyr greater. This could simply reflect 
the correlation of earlier formation (by as much as ~ 1 
Gyr) with high-density cluster environment (e.g. Sheth et 
al. 2006; Rogers et al. 2010). The BCGs also have a flat- 
ter age-luminosity relation, as might be expected if BCGs 
(of all masses) form from dry major mergers of red early- 
type galaxies (with older ages associated with the cluster- 
core environment), which would increase luminosity without 
changing the mean stellar age. In the same way, their CMR 
is flattened. If, on the other hand, the BCGs had formed 
'monolithically', without major merging, they might show 
stronger trends in their stellar populations (age, composi- 
tion) with luminosity/mass, extrapolating the trends seen 
in less massive E/SOs. 

In rest-frame model-magnitude g — r, BCGs are slightly 
redder than non-BCG E/SOs of the same age, whereas in 
spectra-derived g — r, they lie on the same (approximately 
linear) relation of colour to age. We conclude that the 
(~ 0.01 mag) redder model-magnitude g — r colours of the 
BCGs in the CMR and CcrR can be explained simply by 
their greater mean ages combined with their flatter colour 
gradients, which would give them a slightly redder model- 
magnitude colour for a given central-aperture (e.g. spectra- 
derived) colour. 

(viii) We re-examine the correlations of colour gradi- 
ent with galaxy properties, dividing by age (star-formation 
redshift). The positive correlation of colour gradient with ra- 
dius, and the anticorrelation with density and 10 log a + M r 
are strongest in younger galaxies but are seen within all 
age intervals. With increasing age, these correlations flat- 
ten somewhat while the mean colour gradients decrease, as 
might be expected for sequential dry mergers in the Di Mat- 
teo et al. (2009) model. 

The decrease in colour gradient with the a residual 
10 log a + M r is not solely due to the known correlation 
of this quantity with stellar age (Forbes and Ponman 1999; 
Bernardi et al. 2005; Gallazzi et al. 2006), as it is seen within 
each age interval. Rather, it seems the colour gradient must 
have independent anticorrelations with the a to luminosity 
ratio and with stellar age. 

It seems the process of spheroid formation sets up a 
met allicity /colour gradient, with a steepness positively cor- 
related with the mass/luminosity of stars and also with a 
large radius and low mass density for the galaxy. Subse- 
quent mergers flatten the stronger gradients while adding 
low-gradient galaxies at the high-luminosity end of the E/S0 
luminosity function. As dry mergers tend to increase size 
this also weakens the gradient-density relation (Figure 23). 
A comparison with colour gradients in much higher redshift 
E/SOs would help to confirm this. The mean colour gra- 



dient in E/SOs decreases somewhat at Md yn > lO n ' 5 M0 
because most of these higher-mass galaxies are the products 
of elliptical-elliptical (dry) mergers (Tran et al. 2005; Naab, 
Khochfar and Burkert 2006). Yet some E/SOs remain today 
with colour gradients more than twice the observed aver- 
age, like those in the 'monolithic' models of Kobayashi et al. 
(2004). 

The position of BCGs on Figures 21 to 24 is that of 
galaxies which have experienced more age and/or merger- 
related colour gradient suppression than even the Zf orm > 2 
E/SOs. On Fig 18 the BCG colour gradient matches that 
of non-BCGs of mean stellar age > 10 Gyr, although the 
BCG mean stellar age is only 7.5 Gyr. The colour gradients 
in BCGs could be effectively 'super-aged' by several Gyr by 
their environment and their position in the path of infalling 
cluster spheroidals. 
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